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if AS t and At be the corrections to the parallactic coefficient and 
solar parallax 


At = + 


14*2 




14-2 


-37 88 a^ + 1366 AS a ] 


3074 3074 

= —- A/j 4" --ASjy 

n*5 32-0 


It should be mentioned that Newcomb in 1876 attributed 
the apparent increase of the variation to systematic errors of 
observation varying with the age of the Moon. 


On Oscillating Satellites. (Second Paper.) 
By H. 0. Plummer, M.A. 


1. In the earlier paper {Monthly Notices , vol. lxiii. p. 436), to 
which the present is a continuation, periodic orbits which are 
possible in the neighbourhood of centres of libration have been 
investigated with a degree of approximation corresponding to 
the second order of the relative coordinates. The results appear 
to suggest that a still closer approximation may be made with 
some advantage, although the discussion relates only to a small 
part of the restricted problem of three bodies, and has appa¬ 
rently no direct practical bearing. The notation of the former 
paper will be followed generally, but some changes will be 
admitted. The three collinear centres of libration alone will be 
considered here. 

2. The general equations of motion of the satellite are taken 
in the form 


where 


x—2y = O io , y + 2X = 0 01 




O'rj — 


8 r+ *fl 

dx r dy* 


Now at the collinear centres of libration y = o and £2 IO = = o; 

and O rj = o when s is odd. Hence in relative coordinates the 
equations of motion are 

$-2? = 

r) + 2^=zQ l02 T]-{-Q ll ^ r i + i(3& 22 £; 2 0 + Oo 4 ?? 3 )+ • • • 


so long as the satellite remains within the region over which the 
expansions are possible. It is known that periodic solutions of 
these equations exist, and these can be appropriately expressed 
by Fourier series of the form 


£ = a Q + a z cos mt + a 2 cos 2mt-j-a 3 cos ^mt+ . . 

1 7= h x sinm£ + 6 2 sin 2mt-\-b 3 sin 3m£+ . . 
the orbits being symmetrical with respect to the £-axis. 
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If these expressions for 5 and 17 be substituted in (1) and 
powers and products of the circular functions be expressed in 
terms of circular functions of multiples of mt , there will result 
an identity between two cosine series and an identity between 
two sine series. Now it is known that a lt b x are of the first 
order, and a OJ a 2 , b 2 of the second order. The order of the other 
coefficients then corresponds to their suffixes. For if the series 
(2) be supposed to terminate at the argument imt, the series 
obtained by substitution in (1) must terminate at the same 
argument; since it is easy to see that the argument i'mt (i/ > i) 
can only arise in connexion with coefficients of order higher than 

Hence the two identities will give (2^+1) equations involving 
the (2^+1) coefficients and m; and these equations, together 
with the Jacobian integral, will suffice to furnish a solution to 
the order of approximation considered. 

3. The immediate object is to find a solution which shall be 
valid to the third order. The expansions of (1) and (2) may 
therefore be supposed to end with the terms actually written 
down. Then when (2) are substituted in (1) the following 
seven equations are obtained by neglecting terms of higher order 
than the third, and equating the coefficients of the sines and 
cosines of the different multiples of mt: 


O — H 2 o a o H" i^ 3 o a i 2 + 4D I2 6 i 3 
o = (m 2 + £l 20 )a x + 2mf>! +L 
o — (m 2 + Q 02 )b 1 + 2ma I + M 
o = (4m 2 + ^ 2 o)«2 + 4^2+J0 3O a I 2 — ±£l J2 bS 
o — (4m 2 + 0 02 )6 2 + 4ma 2 + |0 I2 a I ^ I 
o = (9m 2 + Q 2 Q )a 3 + 6 mb 3 + A 3 ... 
o = ( 9 ? 7 i 2 4 - 0 02 ) 6 3 + 67 na 3 + B 3 ... 


(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


The quantities L, M, A 3 , 
represent the expressions 


B, are of the third order and 


L = i0 30 a I (2ao + a a ) + iO I aM 2 +-S^4 0 a I 3 f 
M = •|H I2 (2& 0 & i — t 6 2 ) 4 ~iH 22 & I 2 5 I + -|n o4 6 I 3 

A 3 = ■—■Cl 30 ct l ct 2 •5H I2 6 I f> 2 -g-0 22 ^ I ^ I 2 

B 3 := - iH I2 (<x I 6 2 + a 2^i) “b-g ^22 a i 2 b 1 — 


(4') 

(5') 

(S') 

(9') 


4. An eighth equation is required, and this can be obtained 
from the Jacobian integral. Now 

£ = — ma 1 sin mt — 2ma 2 sin 2 mt — 3 ma 3 sin $mt 
i) = mbj cos mt -f 2mb 2 cos 2 mt -f pnb 3 cos $mt 


and 

2CI = 20 0 + a 20 ? 2 + D 02 Ty 2 +|n 30 ? 3 + n i2 ^ 2 

+ t T 2^ 40 ? 4 4- 

I 
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In substituting in the equation 

”V 2 = C 0 —C + 2(0—O c ) 

terms of the fourth order are to be retained, and there results an 
equation involving cosines of the first four multiples of mt. 
Their coefficients must clearly vanish, since the integral must be 
satisfied throughout the motion, and the four relations thus 
obtained must evidently be involved implicitly in the equations 
(3)—(9). There will thus remain an equation involving the 
energy constant which takes this form 

o=C 0 —C +£«i 2 (0 2 o— m 2 ) -f ±bS(fl 02 - m 2 ) + N ... (10) 

where N stands for this expression of the fourth order 

N = ^2 20 (2 a 0 * + a 2 2 )+i£l 02 bs 

+j& 30 aS(2a 0 + a 2 ) aj> z + 2 a t b 2 ) 

+ 3W*i 4 + -h&^bs + 3iA>A 4 “ 2 m*{a* + b 2 ) 
= + 4m 2 (a 2 2 + 6 2 2 )—4ma 2 6 2 ... (10') 


The eight equations (3), (4), . . . (10) are sufficient to determine 
m and the seven coefficients a Q , a I9 . . . 

5. Since (3), (6), and (7) give & OJ a 2 , and b 2 in terms of a ti b z , 
and m, the latter quantities are determined by (4), (5), and (10). 
An approximate solution is obtained by neglecting L, M, and N. 
This will give 

a i _ —2m _ m 2 + Q 02 ./ v 

Y z ra 2 -f-£i 20 —2m 

o = (m 2 + 0 20 )(m 2 + 0 02 )—4m 2 . (12) 

af = (m 2 + 0 02 )a 2 ; b 2 = (m 2 + O ^ 2 ... (13) 

C 0 —C = (?7i4—0 20 0 02 )a 2 .(14) 

Then, if L, M, and H be reinstated, it will be necessary to 
replace a z by a z + la I9 b z by 5 I + 86 I , and m by m + £m, the three 
increments being clearly of the second order. Hence 

o = (m 2 + £l 20 )( : a I -\-2m8b z + 2(ma z -\-b z )8m + Ij 
o = 2m^ I 4-(m 2 4-0 02 )^ I + 2(m5 I + a I )8m + M 
o — (O zo — m 2 )a z 8a z + (£l 02 — m 2 )b z 8b z — m(a x 2 + b z 2 )$m+l$ 


the quantities L, M, N and the coefficients in these equations 
being calculated with the approximate values of a„ b Zi and m. 
Forms more convenient for calculation can be deduced by the 
help of (11), (12), and (13). Thus if 


cbj ^ ^ 

— — c x , — 4 7 

o z a z Oj 
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the following form an equivalent set of equations : 

(ra4— 0 20 Q 02 )Sm = ’L+a " 1 M) . (15) 

4 raSj; = (Q^ —12 02 )8w& + &r x L—o^ -1 M. (16) 

= m 2 (0 20 — ^02)^1—™(2m 2 + 0 20 + 0 02 )8m + a ~ 2 N 

... (17) 

The first of these determines approximately the variation of 
the period from its limiting value; the two latter give suc- 
sively and d 2 , whence Sa x and 8 b T are obtained immediately. 
The coefficients a 3 and 6 3 are found from (8) and (9), and thus 
the solution to the third order is complete. 

6. To the derivatives of the force function given in the 
former paper must be added those of the fourth order. Their 
values are given by 

3^40 = ~ 60 22 = 80 04 = 72 (fir-s -f vr' 5 ) 

A table of numerical values for the case in which up = 10, 
nv == 1, will now be given, many of the numbers being taken 
from the earlier paper : 


(m4—0 20 0 02 )8 2 



Family a . 


Family b . 

Family c . 

ri 

Q7I75I 


1*34700 


0-94693 

r 2 

0*28249 


0*34700 


1-94693 

Co 

3 65292 


3-53418 


317322 

n 2Q 

+ 13*9878 

+ 

6-0955 

+ 

3*1660 

{^02 

- 5*4939 

- 

1-54/8 

- 

00830 

n 30 = — 2 2 

+ 65*076 

— 

39-279 

+ 

6*822 

3^40 — — 6x^22 — 8 jQq 4 

+ 3982*86 

+ I 3 I 5'79 

+ 86*202 

m 

2*6082 


I*6763 


1*0706 

TtlA — £ L 2 qCIq 2 

123*13 


I7*330 


1*5766 

ar z . »i 

- i*I 439 

— 

1*1234 

- 

1*0312 

a~ l . b z 

+ 4*5597 

+ 

2*984I 

+ 

2*0766 

a ~ 2 . a Q 

+ 10*5685 

— 

5*1400 

+ 

0*5887 

ar 2 . a ? 

- 4*1352 

+ 

2*626l 

- 

0-5487 

a ~ 2 . b 2 

— 1*9209 

+ 

1*5800 

— 

0-2895 

a ~3 . L 

+1234*6 

+ 

73*93 

+ 

2-8618 

or 3 . M 

+ 3493'* 

+ 

47*45 

+ 

1-4737 

ar 4 . N 

+ 2981-3 

- 

11874 

- 

2*4l69 

, A 3 

— 2044*8 

- 

288*49 

- 

8-3943 

a ~3 . B 3 

-2190*5 

- 

225*81 

— 

65525 

a* -2 . dm 

“ 153*79 

- 

2*832 

- 

0*03708 

a ~3 . 5 ai 

- 92*47 

— 

6*227 

- 

O367O 

a 3 , 85 j 

+ 224*79 

- 

3*652 

- 

05638 

a ~ 3. a 3 

+ 20*19 

+ 

7109 

+ 

0*4530 

a ~3 . 6 3 

+ 3364 

+ 

6*500 

+ 

o *3559 


1 2 
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7. In the former paper the results for orbits of family a were 
compared with four orbits given by Professor Darwin. The 
comparison can now be repeated with improved figures. The 
values of the 1 energy constant C and the parameter a are 
given by 

11 C= 40*0, 39’o, 38*5, 38*0 

loga= 8*06428, 8*47046, 8*54705, 8*60356 


The general equations for the family are 


£ = + io*5685a 2 —(ri439a + 92*47a3) C os mt—^i^a ' 2 cos 2 mt 

-i-2o*i9a3 cos 3 rat 


7 ) = 


(4*5597^ +224*79^3) sin mt— 1 -9209a 2 sin 2mt 

+ 33*64^3 sin ynt 


m = 2*6082 —153’79^ 2 


which become in the four cases considered 

(a) £ — +*0014 —*0134 cos mt —*0006 cos 2m£-f*oooo cos 3 mt 
71 = +‘0.533 sin mt—' 0003 sin 2m£ + *oooi sin 3 mt 

m = 2*58^ 

(/?) £ = +*0092 — 0362 cos mt—— *0036 cos 2ratf + *ooo5 cos 3 mt 

7) = +*1405 sin mt —*0017 sin 2mt + 'ooog sin $mt 

m = 2*474 

(y) £ = +*0131 —*0443 cos mt — *0051 cos 2mt -{- *0009 cos 3 mt 

r]=: +’1705 sin mt —*0024 sin 2mt + 'ooi$ sin $mt 

m — 2*418 

(r) £ = + *0170—*0519 cos mt — *0067 cos 2m£ + *ooi3 cos $mt 
7 ] = -f- 975 sin mt —*0031 sin 2ra£+*oo2 2 sin $mt 

m — 2*361 

The following table gives the points where the orbits cross 
the axis of x (f Q and £j), the points (£', ?/) for which mt = -Jtt 
and the angle ( 2Trm ~ I ) through which the system turns in one 
period of the satellite. In the comparison numbers rf is the 
ordinate corresponding to the calculated i!. 


mt - 

0 

7r 

JL 

' 



to 

t> 

r f ' 

V 


(“) 

— '0126 

+ -0142 

+ *0020 

+ *0532 

I 39 ° l 

D. 

— •0125 

+ •0147 

»» 

+ -0530 

1383 

0 ) 

-•O3OI 

+ O4I3 

+/OI28 

+ •1396 

1455 

D. 

-•O3O4 

+ •0413 

1 f ‘ 

+ •1395 

1457 
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mt— 

O 

IT 

4^ 




£0 

£1 




(y) 

-*0354 

+ •0514 

+ *0182 

+ *1690 

H &9 

D. 

-0361 

+ *0523 

?» 

+ •1691 

1496 

(8) 

-•0403 

+ -0609 

+ *0237 

+ •1953 

152-5 

D. 

-•041 

+ •060 

»5 

+ •197 

(not deter- 


mined) 


The limiting value of 2nm~ z for an infinitesimal orbit is 
i 38 °*o. The agreement shown in the above table is much better 
than in the corresponding table of the previous paper (vol. lxiii. 
p. 442), and may be considered satisfactory. 

8. For the two orbits of family b given by Professor Darwin 
the values of C and a are given by 

110 = 38-5, 38-0 

log a == 8*64750, 8*83116 


The general equations for the family are : 

i; = —5*1400 a 2 — (i*i234a + 6 , 227a3) CO s mt + 2'626ia 2 cos imt 

4 7*109^3 cos 3 mt 

7} = (2‘9841a—3*652^3) sin mt+i -5800a 2 sin 2 mt 

+ 6*500^3 sin ynt 

m = 1*6763 — 2*832 a 2 


which become in the two cases considered 

(a) £ = — *0101 — ‘0504 cos m^+*oo52 cos 2m£4-*ooo6 cos 3 mt 
rj = +-1322 sin m^+*oo3i sin 2m£4-*ooo6 sin 3 mt 

m = 1*6707 

(fl) f = —*0236—*0781 cos ratf + *oi2i cos 2mtf+*oo22 cos $mt 
r] = 4- '2012 sin mt 4*0073 sin 2mt + 'o&2o sin 3 mt 

m = 1*6633 


The copiparison is made as in the preceding paragraph 


mt- 

O 



JL 



f. 

£1 



27rm 1 

(«> 

-0547 

. +‘0449 

-•0153 

+ •1316 

215-5 

D. 

-•0551 

+ ‘0445 


+ •1312 

1213-9] 

t*> 

- 0874 

+ *0644 

-•0357 

+ *I99I 

216*4 

V. , 

-•0875 

+ *0679 


+ *I987 

208 


The limiting value of 2* m~ l for an infinitesimal orbit is 214°* 8. 
There is therefore a discrepancy between the two methods in 
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regard to the periods of the orbits, Professor Darwin's results 
indicating a decrease of the period as the orbits increase in size. 
This is not supported by my figures. The value [2i3°*9] really 
belongs to a neighbouring orbit which is not quite periodic, so 
that this number is not strictly comparable. On the other hand 
the orbit (/ 3 ) contains in the value of £ x a second discrepancy, 
which is greater than the evidence based on other comparisons 
would lead us to expect. On the whole it seems quite possible 
that the determination of this orbit by the method of quadratures 
has not been so accurate as in other cases. 

9. In Professor Darwin's results the range of the energy con¬ 
stant has the limits given by 11C = 40*0 and 11C = 38*0. Now 
orbits of family c only begin to appear after 110 = 35*0, and 
therefore orbits of this family for comparison are wanting. But, 
although the results cannot be submitted to this kind of check, 
it seems worth while to apply them to drawing a few new orbits 
of this family. The general equations are 

£ = +*5887 a 2 — (ro3i2 ^ + *3670 a 3 ) cos mt— *5487 a 2 cos 2 mt 

+ •4530 cos $mt 

r\ = (2’0766 a— *5638 a, 3 ) sin mt — '2895 a 2 sin 2 mt 

+ *3559 a3 s * n 3 mi 

m = 1*0706 —*03708 a 2 . 

The selected values of 0 and a are given by 

110=34-5 34-0 33-5 33-0 

log a = 9-18439 9-35 8 86 9 ' 4 S 434 9 ’ 5 2 °43 

The equations of the four orbits then become 

(a) £== + *0138—*1590 cos mt —*0128 cos 2mt + *0016 cos $mt 

7 ] = + *3155 s i n m t —*0068 sin 2ra£ + *ooi3 sin 3 mt 

m = 1 0697 ; 27 rm“ I = 336°*55 

(/ 3 ) £ = + *0307 — *2400 cos mt— *0286 cos 2ratf-poo54 cos 3 mt 

rj = +*4678 sin mt —'0151 sin 2mt + *0042 sin 3 mt 

m = 1*0687 > 2nm~ I = 336°*85 

(y) £ = -f *0477 —*3020 cos mt — *0445 cos 2ra£-h*oio5 cos 3 mt 

7 i = +*5782 sin mt— ’0235 sin 2mt + *0082 sin $mt 

m = 1*0676; 2nm~ I = 337°*2o 

(8) £ = + *0647 — *3552 cos mt — '0603 cos 2w£-poi65 cos $mt 

7 j = + *6678 sin mt —*0318 sin 2mt + *0130 sin 3 mt 

m = 1*0665 ; 27 TW 1 = 337°*55 
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The curves corresponding to these equations are represented 
in the accompanying plate (Plate 2), in which the unit of length 
is 5 in., this being the scale on which most of Professor Darwin’s 
orbits are reproduced. 

10. A comparison of §§ 7 and 8 with §§ 10 and 11 of the 
previous paper will show the degree of improvement obtained for 
families a and b by increasing the order of approximation, and 
give some idea of the convergence (in a practical and not tech¬ 
nical sense) of the series. It has now been found possible to 
calculate the variation of the period, an element which, with the 
lower approximation, is supposed constant for each family. In 
every case the period is found to increase with the scale of the 
orbits, the variation being most marked in the case of family a 
and least in the case of family c. As regards the shape and 
position of the orbits, the general features which characterise the 
departure from the small elliptic paths surrounding the centres 
of libration are similar in all cases. The mean position of the 
satellite tends always in the direction of the nearer body of finite 
mass, while at the intersection of the orbit with the line of 
syzygies in the same direction, the curvature first decreases and 
finally changes sign. These features are shown, though not in a 
very marked degree, by the orbits of family c drawn in this 
paper. Where comparison was possible with the results of Pro¬ 
fessor Darwin’s researches it has been made, and in one case 
among the six orbits available a discrepancy is found which is 
pointed out in § 8. The method of the present paper is extremely 
simple and elementary ; yet it might be extended without serious 
difficulty to higher orders of approximation if this were desirable. 
The next order in particular might easily be added. But it 
has not seemed likely that anything of value would be gained 
thereby. A series of orbits of oscillating satellites has been 
found by M. Burrau (A.N. 3251) to be limited by a “trajectory 
of ejection.” The real obstacle in the way of following up to this 
limit the continuation of the class of solutions considered in 
this paper is the limited range over which the force function can 
be expanded in ascending powers of the relative coordinates. 
Whether this difficulty can be surmounted is a question which 
must be left for further consideration. 

University Observatory , Oxford : 

1903 December 7. 
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